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Abstract. The density of zeros of the partition function of the Ising model on a class 
of treelike lattices is studied. An exact closed- form expression for the pertinent critical 
exponents is derived by using a couple of recursion relations which have a singular behavior 
near the Yang-Lee edge. 
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It is well known that Ising model on finitely ramified fractal lattices [1-2] can not 
display a phase transition at any finite temperature. This means that the free energy 
density /(T, H) of such a model provides an analytic function of temperature T (T > 0) 
and of a real magnetic field H. Yang and Lee [3-4] were the first who pointed out that 
the free energy, even of a finite Ising system, can exhibit a singular behavior if one allows 
the field H to become complex H = H' + iH" . They showed that all of the zeros of the 
partition function of a nearest-neighbor ferromagnetic Ising model lie on the unit circle in 
the complex activity plane y = exp{—2H/kgT) = exp(— 2/i). In the thermodynamic limit 
these zeros are expected to condense, yielding to a limiting density of zeros g{h", T). For 
any fixed temperature there exists a pair ±iHQ{T) of zeros lying closest to the real H axis, 
commonly referred to as the Yang-Lee (YL) edge singularities. These limiting nonanalytic 
points exert the most direct influence on the behavior of /(T, H) for real T and H, and it 
is very important, therefore, to understand the nature of the associated singularity in the 
density of zeros near the YL edge. 

The calculation of the limiting density of zeros is generally a highly nontrivial problem 
and little is known about its behavior near the YL edge, let alone about its actual form 
on the whole region of interest. It is widely accepted, however, that the density of zeros 
near the edge exhibits a power-law behavior: (7 ~ \h" — ho{T)\^ , as h" ho(T), where 
a is the YL edge singularity critical exponent. Since in this case exists only a single 
relevant variable [5-8], all other YL critical exponents can be expressed in terms of a. In 
particular, the correlation length critical exponent i^c, which describes the spatial decay of 
the two-spin correlation function near the edge, ^yl ^ \h" — hQ{T)\~'^'' , can be related to a 
{a = dvc — 1, with d being the dimensionality of underlying lattice). It has been suggested 
[6-8] that the value of this exponent depends only on the dimensionality of the lattice and 
is independent of temperature for all T above the critical temperature. In (i = 1, the 
problem can be solved exactly which yields to cr = — 1/2 for T > 0. In spite of the fact 
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that g{h" ^ T) is not known exactly for the two-dimensional Ising ferromagnetic, it has been 
predicted a = m. d = 2 [9-10]. It is believed that a maintains its mean-field value 

(7 = 1/2 above the upper critical dimension cZ = 6 [7]. Only approximative analytical and 
numerical results are available for the intermediate values of space dimensionalities [7-8] . 

The original YL circle theorem [3-4] is independent of the topological structure of the 
lattice and should also apply to the appropriate models of certain non-homogeneous sys- 
tems (diluted ferromagnetics, for example). In order to better understand certain aspects 
of their critical behavior, few studies [11-12] of the density of zeros of the Ising ferro- 
magnetics on a variety of deterministic fractals have been performed so far. It is shown 
that in this case the density of zeros exhibits a scaling form near the YL edge which is 
more complicated than a pure power law. Using the exact recursion relations nontrivial 
values of the edge exponents a were found. Most of these values were estimated only 
numerically, due to the absence of a pertinent fixed point [12] when one uses the standard 
decimation transformation approach. To overcome this difficulty, one can try to construct 
a more liable renormalization-group scheme possessing the appropriate YL fixed point - 
the block-spin transformation, for example. Unfortunately, such an approach is usually 
much more elaborate, and, as soon as we know, only few exact calculations are presented 
so far [11]. In this Letter we develop an approach which allows us to obtain the exact 
closed-form expressions for the edge exponents of Ising models on finitely ramified fractal 
lattices. It turns out that relevant recursion relations themselves have a singular behavior 
near the YL edge fixed point, which is reminiscent of a singular structure of corresponding 
mappings appearing in some recently studied Gaussian models [13-15]. 

Consider the nearest-neighbor ferromagnetic Ising model with a uniform magnetic 
field on a family of treelike lattices shown in Fig. 1. Each lattice of this family can be 
labeled by an integer p, which represents the maximum coordination number of the lattice. 
This parameter can take all values from 2 to oo, with p = 2 corresponding to a simple one- 
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dimensional chain. Let us note that each of these lattices have a finite fractal dimension, 

D = In p/ In 2, which make them very different from a Cayley tree lattice. The partition 

function of the Ising model of an r-th order lattice can be written as a combination of 

only three partial partition functions: z[''^ = (+,+), Z^''^ = Z(^) (-,-), and Z^""^ = 

Z^'^\+, — ), where Z^'^\+, +), for example, denotes the partial partition function with two 

outer Ising spins [see Fig. (lb)] being fixed in the 'up state'. It is convenient to express 

the resulting recursion relations in terms of two reduced variables Z2 = Z{+, —)/Z{+,+) 

and zs = Z{—,—)/Z{-\-,-\-). Thus, for a lattice of index p, we have found 

J _ z l{z2 + z^Y-^ + y^-^zl{l + z^Y-'' 
zl{z2 + z^)P-'^ + yP-\l + zz)P-'^ ' 
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with y = exp(— 2/i). To obtain an r-th order partition function z^^^ {i — 1,2) one has to 
iterate the above recursion relations r times, starting with the following initial conditions: 
-^2°'' = ^iid z^^ = xy, where x = exp(— 2K) (here K = J/k^T > stands for standard 
ferromagnetic interaction strength). It is easy now to express various derivatives of the free 
energy density in terms of these variables and their derivatives with respect to a suitable 
variable. For example, the average magnetization M^*") per spin is given by 

I + Z2 + 2Z3' ^ ^ 

where we have omitted the iteration index r on the right-hand side of (6), and jV^''^ 
represents the number of sites of the underlying lattice (A/'*^^-' = + 1, for a lattice of index 
p). In the above formula ti, t2, and ^3 denote the scaled derivatives t\^^ = {dZ^^^ /dh)/z[^\ 
i = 1, 2, 3, and they also can be calculated by a recursive procedure. In a similar way, one 
can express the two-point correlation function G^'^^ in terms of z^"^ and z^"^ 

GM=<S.S.>-<5,><5.>=^^ii|^. (3) 

The density of zeros g{h" , T) is known [3] to be given by the limiting behavior of the real 
part of M/tt as h' —>■ 0, so that we may focus our attention here on the asymptotic behavior 
of the above variables near the edge. 
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As it has been emphasized, for any T > there exists a strip \H"\ < Hq{T) inside 
of which the density of zeros vanishes. The hmiting value Hq{T) can be determined 
numericaUy, by a study of the magnetization (2) and the correlation function (3). Such 
an analysis reveals that, for H' = and H" — > Ho{T), Z2 and iterate toward the fixed 
point 

2(p-l) p-1 

4 = {-yo) " , 4 = -{-yo) " , (4) 

where i/q = exp(— 2z/io) depends on x (for x = 1/2, we have found ho = 0.328 648 956 . . . 
and ho = 0.152 769 589... for p = 3 and p = 5, respectively). As it can be verified, 
however, recursion relations (1) become singular at this point, leading to a failing of the 
common fixed-point analysis. This puzzle is quite similar to the one encountered recently 
in the studies of critical properties of ideal polymer chains on fractal space [13-15]. In 
particular, both numerical and analytical examination of the above recursion relations 
reveals the existence of an invariant 'line' zs — 2:3(2:2), all points of which are attracted by 
(4). An asymptotic equation of this line, valuable near the fixed point, can be expressed 
as an expansion over the small variable Sz = z^ — z^ (note that at the fixed point one has 
^2 = -sf)- Thus we obtain 



Zz^zl^Cxsjzl-Z2^ C2{zi -Z2)+0[{zi- 23)^/^], (5) 

where the coefficients ci and C2 are given by 



1 + VI + 4p 
ci = ^ , 

(6) 

_ p2 _ 2j9 - 1 + (3j92 - 6p - 1)4 + [p^ -3p-l + (3p3 - 2p2 -7p- l)z^]ci ^ ' 
"""" 2^2^(1 + 23*)(2 + ci + 2pci) ■ 

Using the above formulas one can show that along the invariant line, near the fixed point, 

5z renormalizes according to the law: 5z' = 25 z/ (1 -\-2p -\- ^1 + Ap). This enables us to 

describe the way in which z^ approaches its fixed point value: 



In a similar way we can extract the leading asymptotic behavior of dz^^ jdh and dz^^ /dh 
at the edge. Indeed, it is easy to see that these derivatives satisfy a couple of linear 
recursion relations, an analysis of which reveals the following behavior 



(r) 



dh 



dzo 



(r) 



dh 

h=iho 



~ 2\ (8) 

h=iho 



In the above established asymptotic relations we have supposed h = iho. In fact, 
they also hold for a finite but very small value of 5h = h" — Hq (0 < <^ 1), provided 
^ ^ ^0) where tq ^ 1 is the number of iterations one can make along the invariant line 
before going away from it (starting with the initial conditions in which y = yo exp{—2iSh)) . 
This number depends on the value of Sh and can be estimated from the following obvious 

Sh. Thus, taking into account (7) and (8), we find 



relation: z^ ° {Sh) ^ z^ -\ 

X^o {Sh)~^, where A = l + ^/T+~4p. This enables us to express the YL correlation length 
as a function of 5h: ^yl ^ 2''° ~ exp(— ln2 ln(5/i)/ In A), i.e. 

, In 2 In 2 

independent of the temperature. We have also calculated the critical exponent Vc for 
several values of p, by using a numerical study of the correlation function (3) , and we have 
found an excellent agreement with (9). Let us note further that for p = 2 the value of 
i^c = 1/2) coincides with the known value for the Ising ferromagnetic chain, while for 
p = 3 our exact value of Vc is very close to the one that has been obtained earlier by using 
somewhat different numerical approach [12]. What is perhaps more interesting, our values 
of z^c for p = 2, 3 are the same as those found for the correlation length critical exponent 
z/q of a simple zero-field Gaussian model on the same lattices [13-15]. It is tempting, 
therefore, to see whether this coincidence persists for all values of p. Our calculation [16], 
which follows the lines of Ref. [13], shows that this is, really, the case {uc = ^a)- 

In a similar way we have found: t^"^ ~ t2^^ ~ t^''^ ~ A*^. Taking 

h=iho h=iho h=iho 

into account these results and (2), we can write M^'^°^ ~ t^^°^ /H^^°^ ~ {X/pY°, which 
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yields to 



M^g^Sh", a=— , - 1. (10) 



It is interesting to note that this expression can be rewritten in the form of an expected 
scahng relation: a = Duc — 1. As it can be seen from the above formula, a — (j{p) is an 
increasing function of p and its limiting value, a{p oo) = 1, exceeds the above mentioned 
mean-field value a = 1/2 for the homogeneous lattices. Note that a{p) changes its sign 
at p = 6, which means that singular part of the density of zeros diverges near the edge if 
p < 6, whereas it vanishes for all p > 6. As can be noticed in Fig. 2, where we plotted the 
density of zeros as a function of h" for several values of the lattice parameter p, the zeros 
are not homogeneously distributed on the unit circle but instead have a striking Cantor-set 
structure with many gaps, especially for the low values of p (except, of course, for p = 2). 
As p increases, this fractal pattern become less and less pronounced, and it seems that, for 
any small but finite h' , it acquires a rather regular limiting form (see Fig. 2. for p ^ 1). 
For a fixed and finite temperature, the value of the YL edge decreases with p and one can 
show that ho ~ 1/p when p ^ 1, corresponding to the absence of any gap in the density 
of zeros, in the limit when the lattice coordination number tends to infinity. On the other 
hand, for a finite p the density of zeros exhibits a scaling form near the edge, similar to 
the one that has been discussed in details in Ref. [11-12]. 

In conclusion, we have presented exact results for the density of zeros in the complex 
activity plane for a family of treelike lattices. We have found that our recursion relations 
have a singular behavior near the edge. The technique we have used here to extract the 
leading singularity can be, in principle, extended to all other finitely ramified fractal lattices 
and to other classical spin models. 
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Figure Captions 

Fig 1. (a) First two stages in the iterative construction of a p = 4 treelike fractal 
lattice, (b) Schematic representation of an r-th stage p — 3 fractal lattice. In order to 
obtain a partition function zl^^\Si, S2), i = 1,2,3, one has to fix the states of any two 
outer Ising spins Si and S2 (open circles) and perform summation over the states of all 
remaining spins [including the central spin (black circle)] 

Fig 2. Density of zeros g{h" , T) as a function of the imaginary field h" for several 
values of p. In all cases temperature corresponds to a; = 1/2 and the real part of the field 
is h' = 0.005. 
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